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Abstract: The shear stress generated by the wind on the land surface is the driving force that results in the 
wind erosion of the soil. It is an independent factor influencing soil wind erosion. The factors related to 
wind erosivity, known as submodels, mainly include the weather factor (WF) in revised wind erosion 
equation (RWEQ), the erosion submodel (ES) in wind erosion prediction system (WEPS), as well as the 
drift potential (DP) in wind energy environmental assessment. However, the essential factors of WF and 
ES contain wind, soil characteristics and surface coverings, which therefore results in the interdependence 
between WF or ES and other factors (e.g., soil erodible factor) in soil erosion models. Considering that DP 
is a relative indicator of the wind energy environment and does not have the value of expressing wind to 
induce shear stress on the surface. Therefore, a new factor is needed to express accurately wind erosivity. 
Based on the theoretical basis that the soil loss by wind erosion (Q) is proportional to the shear stress of 
the wind on the soil surface, a new model of wind driving force (WDF) was established, which expresses 
the potential capacity of wind to drive soil mass in per unit area and a period of time. Through the 
calculations in the typical area, the WDF, WF and DP are compared and analyzed from the theoretical basis, 
construction goal, problem-solving ability and typical area application; the spatial distribution of soil wind 
erosion intensity was concurrently compared with the spatial distributions of the WDF, WF and DP values 
in the typical area. The results indicate that the WDF is better to reflect the potential capacity of wind 
erosivity than WF and DP, and that the WDF model is a good model with universal applicability and can 
be logically incorporated into the soil wind erosion models. 
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1 Introduction 


The wind is the driving force of soil wind erosion. When the shear stress exerted on the soil surface 
by wind exceeds the ability of soil to resist detachment, soil wind erosion occurs (Skidmore, 1986). 
Therefore, wind velocity is regarded as an important variable in soil wind erosion models 
(Woodruff and Siddoway, 1965; Hagen et al., 1996; Fryrear et al., 1998). Since the shear stress of 
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the wind on the soil surface is an important factor determining the soil loss by wind erosion, the 
relationship between the soil loss by wind erosion (Q) and the shear stress (z) of the wind on the 
soil surface is generally written as Qa(r'?-z,!")r (Lettau and Lettau, 1978), or Qo(t—1;)r!” 
(Sauermann et al., 2001), or Q(z) (Anderson and Haff, 1988; Andreotti, 2004; Duran and 
Herrmann, 2006; Martin and Kok, 2017), where t=pgu~, t=pgu, T and t; expressed in N/m?. Here, 
7: is the critical shear stress when wind erosion occurs, pg is the density of air (kg/m°), and ux and 
u» are the friction velocity and the critical friction velocity (m/s), respectively. When the soil wind 
erosion reaches steady state, Q represents the soil particle flux in per unit width and unit time. 
Previous studies have shown that the number of soil particles (N) that is driven by the wind in the 
unit area and unit time is also positively correlated with (t-t) (Cheng et al., 2007). 

Because of many factors affecting soil wind erosion, previous studies often took the factors, such 
as surface crop residue coverage, crop upright biomass, surface soil moisture, precipitation and 
evaporation, into account when describing the wind erosivity, and called wind erosivity as the 
weather factor (WF) (FAO et al., 1979; Hagen et al., 1996), or the erosion submodel (ES) (Hagen 
et al., 1996). The problem is that WF and ES do not exclusively express wind erosivity, but rather 
combine wind erosivity with the intrinsic properties of the soil (e.g., soil moisture) and surface 
attachments (e.g., crop residues) (Zou et al., 2014). In fact, wind erosivity depends on the wind 
velocity determined by the weather process and does not change with other factors influencing soil 
wind erosion. Roughness elements on the surface (such as plants and ridges) only re-transform the 
wind velocity of the near-surface layer. From the perspective of wind erosion dynamics, it is 
necessary to study the energy relationship between wind and soil, and the factors affecting soil 
wind erosion can be divided into air, surface and soil (Chepil, 1945). During the soil wind erosion 
process, the properties and mechanical attributes of each influencing factor are different. According 
to the attributes that are favorable or unfavorable to soil wind erosion, and the principles of similar 
mechanical properties exhibited during soil wind erosion process, the three categories of factors, 
which are air, surface roughness and soil, represent wind erosivity factors, roughness interference 
factors and soil anti-erodibility factors, respectively (Zou et al., 2015). Corresponding to the three 
categories of factors are wind driving force (WDF), roughness interference force (RIF) and soil 
anti-erodibility force (SAF), respectively (Zou et al., 2014). Studying the energy relationship 
between wind and soil is of great significance during the process of soil wind erosion (Chepil, 
1945). The purpose of establishing the WDF model is to distinguish the wind driving force from 
soil characteristics and roughness factors affecting soil wind erosion. 

Wind erosivity factors include the wind velocity, wind direction, density of air, viscosity of air 
and turbulence (Zou et al., 2015). In the soil wind erosion models, the WDF should be an 
independent variable or a submodel, which is spatially located above the soil surface, and is in 
contact with the surface roughness elements and soil surface. In terms of mechanical properties, 
the WDF reflects the potential for wind to produce shear stress on the soil surface and causes the 
displacement of soil particles away from their original positions (Zou et al., 2014). For a given area, 
soil wind erosion also has a time scale, which means that the WDF must match the unit area and 
time scale. Therefore, when establishing the WDF model, the focus of attention must be the 
accumulative shear stress generated by the wind on the surface per unit area in unit time or a given 
period. In general, the WDF is defined as the accumulative shear stress generated by wind on a 
surface per unit area in one year (Pa h/a). 


2 Theoretical foundation and model establishment 


The information available in previous studies on the WDF was fragmentary and insufficient to 
support the establishment of the WDF model (Zou et al., 2014, 2015; Zhang et al., 2018). It was, 
therefore, necessary to start from the basic assumptions and their verifications. 


2.1 Hypothesis and verification 


2.1.1 Hypothesis 1 
In each wind erosion event, the wind velocity always increases from low to high, gradually reaches 
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the maximum, and then reduces. Also, the wind velocity can be considered to consist of several 
grades of wind velocity. Since u» exhibits different intensities on soil wind erosion, the resulting 
values of t will therefore differ. Consequently, in each wind erosion event, the value of u» must be 
divided into several grades. Previous studies showed that the wind velocity variations in soil wind 
erosion areas in China and abroad are consistent with this assumption (Hagen et al., 1996; Fryrear 
et al., 1998; Yuan et al., 2018), even globally (Dai and Deser, 1999). 

2.1.2 Hypothesis 2 

The WDF model considers that the relationship between the soil loss by Q and z of the wind on the 
soil surface as Qoc(z—1). The model also assumes that when soil wind erosion is at steady state, the 
soil wind erosion rate (kg/(m2.s)) refers to the ratio of the aeolian sand flux (kg/(m-s)) by the length 
of the plot (m), which represents the soil loss by wind erosion per unit area per unit time. In order 
to verify the objectivity of this hypothesis, it is necessary to start by clarifying the following two 
basic concepts. Firstly, Bagnold (1941) initially proposed the aeolian sand flux Qor?*. To avoid 
meaningless result calculated from Qoar”? when r<z;, the forms of Qæ(t"?—tt)t and Qoc(t-1))r!” 
were later proposed (Lettau and Lettau, 1978; Sauermann et al., 2001), these equations are overall 
expressed as Qocr*”. Secondly, In order to reveal in essence the relationship between the shear 
stress of the wind acting on the soil surface and the Q, the Qo(z—1;) form was proposed (Anderson 
and Haff, 1988; Andreotti, 2004; Durán and Herrmann, 2006; Martin and Kok, 2017), where (zz) 
was the effective instantaneous shear stress of the wind acting on the soil surface (N/m?); Q was 
the particle mass of the soil surface that is the instantaneously driven by the wind in the unit area 
under the action of (t-t). When the soil wind erosion is at steady state, the particle mass on the soil 
surface that is driven by the wind per unit area per unit time is completely transported by the wind. 
Therefore, the Q in the aeolian sand flux equations Qor??, Qoc(r!/7—1;""*)t and Qoc(t—-1)r/? divided 
by the length of the plot represents essentially the soil wind erosion rate. 

2.1.3 Hypothesis 3 

In the ‘appropriate time scale’, the soil wind erosion process is considered to be at steady state. The 
‘appropriate time scale' here is a variable, and its specific value is determined by the duration of 
monitoring soil wind erosion. For example, it can be 5 or 10 min when monitoring soil wind erosion 
in a day or a month in a limited area, and it can be an hour when monitoring inter-annual soil wind 
erosion in a large area. The point of this hypothesis is that since the wind is always in a turbulent 
pulsating state by nature, which may induce the soil wind erosion rate to have strong pulsating and 
intermittent characteristics, it is very difficult to strictly define whether or not the soil wind erosion 
is at steady state (Zou et al., 2019). Hence, both the wind velocity is in turbulent conditions, and 
aeolian sand flux have time-scale effects (Duran et al., 2011; Martin et al., 2013; Mayaud et al., 
2017). This means that for a specific turbulent wind, there is a good correlation between wind 
velocity and aeolian sand flux in the average value of a certain time period (Butterfield, 1991, 1998, 
1999; Wang and Zheng, 2013). Therefore, the so-called steady state is actually only limited to a 
relatively short period of time in which the soil wind erosion rate fluctuates within a small range. 
For the study of soil wind erosion process, the length of the time period is the decisive factor in 
determining the ‘appropriate time scale’. Currently, the ‘appropriate time scale’ is usually ranges 
from 5 to 10 min (Udo, 2016; Shang et al., 2017) in a day or a few days, and an hour in the inter- 
annual time period (The Office of the First National Water Resources Census Leading Group of the 
State Council, 2010). 

2.1.4 Hypothesis 4 

Within a given period of time, the accumulative value of (t-t) or z is suitable for the calculation of 
soil loss by wind erosion. This assumption is based on the relationship Qoc(z—z;), which means that 
Q is linear with (t-t) or t. The value of z is only related to the potential capacity of the wind-driven 
soil surface particle mass, which reflects the potential capacity of u» for soil wind erosion. However, 
the occurrence of soil wind erosion depends on whether the value of z (i.e., u» value) is greater than 
Tte (i.e., u» value). When t>7; (i.e., u+>u»), (tt) represents the effective shear stress of the wind 
acting on the soil surface. Therefore, within a given period of time, the relationship between the 
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cumulative value of (z-z;) or t and the cumulative value of soil loss by wind erosion is not only 
linear, but also one-to-one correspondence. Considering the variability of u» with the physical and 
chemical properties of soils, a fixed u» value cannot be given. Also, the WDF is defined as the 
potential capacity for the wind to produce shear stress on the soil surface and to cause the 
displacement of soil particles from their original position (Zou et al., 2014). Therefore, for the 
calculation of the WDF in this study, the accumulative value of t was considered. 


2.2 Principles of the model establishment 


2.2.1 Principle 1 


The WDF model is established for flat surfaces without roughness elements, as constrained by the 
definition of the WDF given above. From the perspective of wind erosion dynamics, the only forces 
acting on soil wind erosion are the WDF, RIF and SAF (Chepil, 1945; Zou et al., 2014, 2015). The 
WDF is an independent driving force of soil wind erosion. RIF is the re-transformation of the WDF 
by various roughness elements attached to the surface, which therefore weakens the WDF and acts 
in opposition with the space located above the soil surface. SAF is a force that the soil itself 
possesses against the WDF and is determined by factors including the soil physical and chemical 
properties, as well as the plant roots. Therefore, when establishing the WDF model, the value of 
aerodynamic roughness length (zo) was considered to be 0.02 mm, in agreement with previous 
studies on a flat surface (Bagnold, 1941; Dong et al., 2001; Liu et al., 2003). Using u» and a fixed 
Zo value when applying the WDF model, the latter is therefore not subjected to the limitations of 
wind velocity observation heights and areas. 

2.2.2 Principle 2 


The WDF model has temporal meaning and can represent the magnitude and duration of the force 
of the wind acting on the soil surface. Soil wind erosion is not only related to wind velocity but 
also related to the action time of the wind. Under similar wind velocity conditions, the soil loss by 
wind erosion increases with increasing wind action time. For similar wind action time, the soil loss 
by wind erosion rapidly increases with the increasing wind velocity. Obviously, the wind velocity 
and wind action time will continuously vary in a certain period of time (hour, day, month and year). 
When calculating the soil loss by wind erosion, the shear stress generated by the wind on the soil 
surface and the wind action time is therefore used in order for the wind velocity and the wind action 
time to be considered simultaneously. 

2.2.3 Principle 3 

The WDF model has a theoretical basis for mechanics and contains easy-to-obtain indices. The 
ultimate goal of establishing the WDF model is to establish a soil wind erosion model that has a 
theoretical basis for wind erosion dynamics. The wind velocity and direction, temperature, air 
pressure and altitude included in the WDF model are the regular measurement constituents of the 
weather station. The density of air can be calculated based on some regular measurement 
constituents including the altitude of the weather station, the air temperature and the air pressure. 
The computation results of the WDF model represent the potential erosivity of the wind on a flat 
surface, considered as the maximum erosivity. When roughness elements are present on the surface, 
the shear stress of the wind shared by the roughness elements is calculated by the shear stress 
partition model. The maximum shear stress, the average shear stress and the shear stress distribution 
actually received on the soil surface can also be calculated by the shear stress partition model 
(Raupach, 1992; Raupach et al., 1993; Raupach et al., 2006; Walter et al., 2012), while these three 
shear stresses on the soil surface are not considered in this paper. 


2.3 Model 


Based on the above assumptions and principles, in a soil wind erosion event, for a certain u»; level, 
the shear stress generated per unit area of the soil surface is: 


Bele (1) 


where u»; represents the i" friction velocity level (i=1, 2, ..., n, i.e., each soil wind erosion event 
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contains n friction velocity levels). Generally, u+; is graded at intervals of 0.034 m/s, which 
corresponds to a level of 1.0 m/s wind velocity interval at the height of 10 m of a weather station. 
When using the wind velocity data observed by the weather station, u»; is calculated using the 
formula u; / us; = ln(z /Zzo)/k> where u; is the wind velocity observed by the weather station; z is 


the height at which the wind velocity is observed; k is the Von Karman constant, which generally 
takes a value of 0.4; and zo=0.02 mm. The mass of soil particles driven by wind per unit area of the 
soil surface is giccu-it; due to the u»; level, where t; is the duration of the ux level. For a soil wind 
erosion event, the total amount of soil particle mass q (kg/m?) driven by wind in the unit area of 


n 
the soil surface is the sum of all qi values such that q = Ba : 
i=1 
In classical mechanics, an impulse represents an increment of the momentum of an object 
subjected to the influence of resultant of external forces. In other words, it is the physical quantity 
that represents the accumulative effect of the resultant of external forces on the object over a certain 
time period, and is the reason for change of the mechanical motion state of the object. Since the 
physical meaning of impulse is completely consistent with the physical meaning that the WDF is 
required to express according to the principles for the model establishment, the impulse can, 
therefore, be used to express the WDF. Theoretically, during the j™ soil wind erosion event, the 
function that the u» changes with time can be given as u-={(T), and consequently 


t=T 

WDF,= f Pg f°(T)aT > where T (s) is the duration of the j" soil wind erosion event. It nonetheless 
t=0 

remains difficult to give an accurate expression of f(T) in practical applications, which renders more 

convenient the application of discretized u». Thus, for a specific u»; level, the WDF; generated on 

the soil surface per unit area is WDF =zit;. In a soil wind erosion event, the total WDF; generated 


n 

on the soil surface per unit area by all levels of u» is I, =) Tit; » Where j represents a given j® soil 
i=1 

wind erosion event (j=1, 2, ..., m) occurring within a certain time period. Hence, if m wind erosion 


events occur within a given time period, WDE-5" lla 
j= 
In practical applications, the value of u» can be accurately determined for a given study area. 
When relating the WDF to the calculation of soil loss by wind erosion, the condition užu» should 
be satisfied for each value of u,; in order to eliminate the impacts of the data of u»<u» on the 
WDF calculation results in the process of calculating the value of the WDF. 


3 Application of the WDF model 


3.1 Application area 


The application area (32°07'12"—49°20'24'N, 75°08'24"-126°34'48"E) of the WDF model is about 
3.95x10° km?. Soil wind erosion occurs the most frequently in this part of China (Ministry of Water 
Resources of China and National Bureau of Statistics of China, 2013). The annual average wind 
velocity in this area is 2.6 m/s. The annual average wind velocity has obviously spatial differences, 
the maximum value of annual average wind velocity is 8.6 m/s (Fig. 1). This area presents an annual 
mean air temperature of about 6.0°C, a mean annual precipitation of 291 mm and a mean annual 
surface snow cover days of around 45 d. 


3.2 Data 


The data used in this paper included the altitude of the weather stations, wind velocity, wind 
direction, air temperature, atmospheric relative humidity, precipitation, ground temperature (0 cm) 
and snow cover. The time series ranges from 1980 to 2016. The wind velocity, wind direction, daily 
air temperature, daily atmospheric relative humidity, daily precipitation and daily ground 


JOURNAL OF ARID LAND 


temperature were derived from 157 weather stations in the application area (Fig. 1). 
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Fig. 1 Application area of the WDF model and the spatial distributions of meteorological stations and annual 
average wind velocity 


There are two sources of wind velocity and wind direction data. One source is the first National 
Water Conservancy Census on Wind Erosion Census group (Zou et al., 2012), with time series from 
1 January 1980 to 31 December 2009. Another source is the National Meteorological Information 
Center, with time series from 1 January 2010 to 31 December 2016. Both data sources provided 
hourly data of u» and wind direction. The accuracy of the wind velocity measurement was 0.1 m/s, 
and the wind direction was 16 directions. 

The data of daily air temperature, atmospheric relative humidity, precipitation and ground 
temperature were from the Dataset of Daily Ground Climate Data in China (v3.0), which was 
provided by the National Meteorological Information Center. The time series ranges from 1 January 
1980 to 31 December 2016. The measurement accuracies of the daily average air temperature and 
the ground temperature, daily average atmospheric relative humidity and daily precipitation were 
0.1°C, 1% and 0.1 mm, respectively. 

There are two sources of snow cover data. The first source was from the long-term snow depth 
dataset of China (1979-2014) (http://westdc.westgis.ac.cn/heihe) published by the Cold and Arid 
Regions Sciences Data Center at Lanzhou. The second was from the accumulative days of monthly 
snow cover provided by the National Meteorological Information Center. These two sources 
provided daily snow cover data. 

In this paper, each year represented a given time period when applying the FDF model. The time 
interval between wind velocity level and wind direction division was set to 1 h. The divided time 
intervals of temperature, atmospheric relative humidity, precipitation and ground temperature were 
set to 1 d. The hourly value of the WDF was calculated by using hourly wind velocity data, and the 
daily value of the WDF was after that derived by accumulating the WDF over 24 h per day. The 
daily actual air density (pg, kg/m*) can be calculated through Equations 2—4. 


Pg = PoloP / (RT), (2) 
P =101.08e 1/7924) | (3) 
T =273+t, (4) 


where po is the air density under standard conditions and is equal to 1.293 kg/m; To is the air 
thermodynamic temperature under standard state and is equal to 273 K; P is the actual atmospheric 
pressure (kPa), which varies with the altitude (Wang, 1984); Po is the standard atmospheric pressure 
and its value is 101.3 kPa; H is the altitude (m); and t is the daily average temperature (°C). The 
data above were all spatial interpolation with Kriging interpolation method. 


3.3 Application results 
In the application area of the WDF model, the WDF generally indicated a downward trend. The 
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Mann-Kendall trend test (Mann, 1945; Kendall and Gibbons, 1990) results showed that 1988 was 
the mutation time point (confidence interval of a<0.05). During the period 1980-1988, the WDF 
rapidly decreased from 234.54 to 179.04 N h/m’. During the period 1989-2016, there was a 
fluctuating slow downward trend. The range of the WDF varied from 103.70 to 154.86 N h/m? (Fig. 
2). 
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Fig. 2 Annual variation of the WDF values and the Mann-Kendall trend test (in the upper-right corner) in the 
application area during 1980-2016 


The WDF in every season showed a downward trend during the period from 1980 to 2016, but 
the WDF value changed significantly in different seasons (Fig. 3). The WDF value in spring was 
significantly greater than in the other seasons, and its decreasing rate was also significantly greater 
than that in the other seasons. The WDF value in spring decreased rapidly from 112.10 to 81.39 N 
h/m? during 1980-1988, and the annual average decline rate reached 3.16 N h/m?. However, the 
annual average decline rates of the WDF values during this period were only 1.25, 1.67 and 0.72 
N h/m? in summer, autumn and winter, respectively. During 1989-2016, the change of the WDF 
value in each season slowed down. The change of the WDF value in spring ranged from 48.85 to 
78.03 N h/m? with an average of 64.26 N h/m?. The average values of the WDF in summer, autumn 
and winter were 31.65, 30.17 and 30.20 N h/m?, respectively. 
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Fig. 3 Seasonal variation of the WDF values in the application area during 1980-2016 


In the model application area, the multi-year average of the WDF was 149.43 N h/m?, the median 
value was 109.00 N h/m?, and about 2/3 of the application area was smaller than the median value 
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(Fig. 4). In the spatial distribution, the high-value areas of the multi-year average WDF value were 
mainly distributed in the north and southwest of the central part of the model application area, 
where the maximum value reached 744.71 N h/m?. Also, the low-value areas of the WDF were 
mainly distributed in southwestern and the southeastern part of the model application area, and the 
lowest value was distributed in southwestern part of the model application area, and was only 4.49 
N h/m?. Based on the Mann-Kendall nonparametric statistical method (Mann, 1945; Kendall and 
Gibbons, 1990), the trend of the WDF was tested. The positive value of the parameter Zc means 
that the general trend is increasing, and the negative value of Zc means that the general trend is 
decreasing. The results indicated that the Zc high-value areas were mainly distributed in the north 
of the central part of the model application area, as well as southwest of the central part of the 
model application area, which meant that the WDF exhibited an increasing trend in these areas. 
Similarly, the Zc negative values demonstrated that the WDF presented a decreasing trend in 
southwestern part of the model application area. 
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Fig. 4 Spatial distribution of annual mean WDF value (a) and Mann-Kendall trend test (b) in the application area 
during 1980-2016 


The WDF values of each season were similar to the annual mean WDF values in spatial 
distribution (Fig. 5). The highest and average WDF values of spring, summer, autumn, and winter 
were 259.55 and 63.41, 166.03 and 27.52, 163.71 and 28.55, and 237.99 and 29.88 N h/m?, 
respectively. 
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Fig.5 Spatial distribution of seasonal mean WDF value in the application area in spring (a), summer (b), autumn 
(c) and winter (d) during 1980-2016 
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4 Discussion 


4.1 Rationality of the WDF model 


There are a large number of factors affecting wind erosion that noticeably increases the complexity 
of soil wind erosion research. Amongst these influencing factors, some are favourable and will 
promote soil wind erosion, while others are unfavourable and will have a rather inhibitory effect. 
The wind is the only factor that plays a positive role in promoting soil wind erosion (Skidmore, 
1986). The WDF proposed herein has a clear physical meaning and contains not only the force of 
the wind acting on the soil surface but also the duration of this force, which is essentially the 
impulse generated by the wind on the soil surface per unit area and its duration. However, the 
physical meaning of other parameters related to the WDF currently in use is not clear. For example, 
WF in the revised wind erosion equation (RWEQ) is expressed in kg/m (Fryrear et al., 1998), and 
the unit of the drift potential (DP) expressing the wind energy environment is in the cubic of per 
nautical mile per hour, that is (n mile/h)? (Fryberger and Dean, 1979). 

The simple use of the wind velocity or u=» to express the erosion ability of wind on the soil surface 
is not comprehensive, because pg changes significantly due to the impacts of air temperature and 
altitude, which may affect the calculation result of t. From 1980 to 2016, in soil wind erosion areas 
in northern China, the monthly average py value ranged between 1.01 and 1.14 kg/m’, the daily py 
value ranged between 0.71 and 1.36 kg/m?. It can be seen that the variation range of py value is 
very large at different time. The WDF model can justly reflect the change of pg and u», as well as 
wind erosivity, but WF and DP cannot. 

In addition, from the perspective of wind erosion dynamics, wind erosivity should be an 
independent variable in the soil wind erosion model, it is only related to the relevant characteristics 
of air (Chepil, 1945; Zou et al., 2014; Zou et al., 2015), but the wind cannot be mixed with certain 
soil characteristics. The purpose of establishing the WDF model is to achieve this goal, and 
overcome the shortage that the wind velocity is mixed with certain soil properties in WF and ES 
(Hagen et al., 1996; Fryrear et al., 1998). 


4.2 Other parameters related to the WDF 


The parameters directly related to the WDF are WF and DP. Among them, the calculation equation 
of WF is (Fryrear et al., 1998): 


WE = Wf x SW xSDx pgg, (5) 
_Na ` 2 
eae as) l (6) 


where Wf is called the wind factor (m/s)?; SW and SD are respectively called the soil wetness factor 
and the snow cover factor which are dimensionless, with values varying within the range of 0-1; 
SW can be estimated based on potential evaporation, precipitation, irrigation, and the number of 
days of precipitation and irrigation (Fryrear et al., 1998), or based on average air temperature and 
relative humidity (Cheng and Cheng, 1980; Dong and Kang, 1994); pg is the dry air density, which 
is usually equal to 1.29 kg/m?; and g is the gravitational acceleration (9.8 m/s”). When calculating 
Wf, U; is the it observed wind velocity (m/s) at the height of 2 m above the ground. When U;<5 
m/s, Wf=0; U; is the critical wind velocity of soil wind erosion (m/s) at the height of 2 m above the 
ground with a value of 5 m/s; N is the number of wind velocity observations in a time period of 1— 
15 d (N=500); Na is the number of days in the observation period (Fryrear et al., 1998, 2000). It can 
be seen from Equations 5 and 6, the core idea of WF is an expression of wind erosivity established 
on the basis of the equation of aeolian sand flux proposed by Bagnold (1941) and Zingg (1953), 
but mixes with the factors other than wind, which caused WF to not really express potential 
erosivity of wind. 

The purpose of proposing DP was to use the reliable ground wind information to study the 
formation and evolution of sand dunes (Fryberger and Dean, 1979), and the theoretical basis of DP 
is based on the equation of aeolian sand flux proposed by Lettau and Lettau (1978). The calculation 
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equation of DP is: 
n 
DP = $ (WV? (Vi -V, It), (7) 
j=l 


where n is the total number of different wind velocity levels, which are divided into five levels; i 
is one of the five wind velocity levels (i=1, 2, +++, 5); Wiis the weight of the i wind velocity level; 
Vi is the i-level wind velocity representative value (n mile/h) at the height of 10 m above the ground; 
V; is the critical wind velocity (n mile/h) at the height of 10 m above the ground; t; is the frequency 
(%) of the i wind velocity level. When calculating DP, different weights are assigned to different 
V; (Fryberger and Dean, 1979; Pearce and Walker, 2005), which simplifies the equation of aeolian 
sand flux proposed by Lettau and Lettau (1978). DP is a substitution index of the wind energy 
environment. It can only indicate the relative strength of the aeolian sand flux and does not have 
the value of expressing potential erosivity of wind. It cannot be applied in the soil wind erosion 
model. 


4.3 Comparison of wind erosion intensity and the WDF in a typical area 


To further compare the differences in the WDF, WF and DP in the expression of wind erosivity, 
WE and DP were additionally calculated in the Taklimakan Desert and its surrounding areas in the 
west part of the soil wind erosion area in northern China. The reason why this area was selected is 
that the terrain is relatively flat, the soil moisture content is generally low, and there are very little 
roughness elements on the ground (e.g., vegetation). In other areas, the soil wind erosion intensity 
can be significantly affected by mountains, vegetation and precipitation, which is not conducive to 
distinguish the impacts of wind erosivity on soil wind erosion intensity. Figure 6 indicates that the 
WDF is the best fitted with soil wind erosion intensity in spatial distribution, DP is moderately, and 
WF is the worst. The west part of the Taklimakan Desert is affected by the loose sediments of the 
dry lake basin and represents a higher level of soil wind erosion intensity, which is slightly different 
from the performance of the WDF. In the calculation of DP, the weight (Wi) of the wind velocity 
level is artificially added, the higher the wind velocity level, the greater the weight, which results 
in the differences between the maximum and minimum values of DP to be improperly enlarged. 
Moreover, the spatial distributions of DP value and soil wind erosion intensity showed opposite 
trends in the most areas of surrounding Taklimakan Desert. 
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Fig.6 Spatial distribution comparison of soil wind erosion intensity with the WDF, WF (weather factor) and DP 
(drift potential) in the typical area (Taklimakan Desert and its surrounding areas). (a), spatial distribution of soil 
wind erosion intensity modified from the literature (Niu, 2017); (b), a scale enlarged view of the typical area; (c), 
(d) and (e), the spatial distributions of the WDF, WF and DP in the typical area, respectively. 
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5 Conclusions 


In the existing soil wind erosion models, the factor or called submodel expressing wind erosivity 
is not independent variable, which makes the theoretical basis of these models unreliable. Each soil 
wind erosion event is a random event that occurs within a certain period of time. The WDF model 
is a variable in the soil wind erosion model, which is a type of random event. The independence of 
random events is a unique and very important concept in probability theory. The theoretical basis 
for establishing any model is that several variables should be independent of each other. According 
to the definition, the WDF does not depend on any other influencing factors of soil wind erosion; 
that is, the WDF model is an independent variable. Therefore, it meets the modeling requirements 
of the soil wind erosion model, and overcomes the limitations that WF and ES present in combining 
wind and soil moisture as well as crop residues (Hagen et al., 1996; Fryrear et al., 1998). 

From the perspective of the earth surface material cycle process, which contains three links of 
erosion, transportation and deposition, erosion is the first link. As a type of soil erosion, soil wind 
erosion should not include the process that the eroded soil particles are transported and deposited 
by the wind. Therefore, it is defined as 'the process of wind force causing the topsoil particles to 
detach the original spatial position’ (Zou et al., 2014). The difference between soil wind erosion 
and soil loss by wind is that the latter is a concept with time and spatial scale (Zou et al., 2015), 
which not only reflects the results of wind erosion but also shows the transportation and settlement 
results of the eroded particles (Zhang et al., 2018). In this paper, based on the theoretical basis of 
Qarz, the WDF model precisely expresses the potential capacity of soil particle mass per unit area 
to be driven by the wind in a certain period of time. 

The WDF model was established under the theoretical framework of the dynamic models of soil 
wind erosion (Zou et al., 2015) and had a clear physical meaning rather than a statistical model 
based on experimental data. Because the WDF, RIF and SAF were strictly distinguished, all the 
three are independent of each other. The WDF can be logically incorporated into the soil wind 
erosion model and is more applicable than WF, ES and DP. The comparison results of spatial 
distributions of the WDF, WF and DP with that of soil wind erosion intensity in the typical area are 
strong evidences (Fig. 6). In practical applications, for an area where u» can be determined, in order 
to eliminate the impacts of the data of u»i<u»: on the WDF calculation results, u» values of =u are 
selected to participate in the calculation of the WDF. The WDF model lays a reliable theoretical 
foundation for the establishment of the dynamic models of soil wind erosion. 
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